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1. INTRODUCTION 

Blazars form up relatively small (~ 1%) sub-group of active galactic nuclei (AGNs), characterised with 
a relativistic jet, viewed at small angle by observer on Earth [1], There are two components in their 
broadband spectral energy distribution (SED): low-frequency one, with the maximum in NIR/optical/X- 
ray band, and high-frequency one, with the peak in gamma-rays. Those two features generally explained 
in terms of synchrotron and inverse Compton mechanisms of emission [2]; corresponding peak frequencies 
indicated with Vp eak and v^ ak . Blazars divided into high/low-frequency peaked, according to the z+ eafc : high 
synchrotron peaked (HSPs) ones with the z+ eafe > 10 165 Hz, and low synchrotron peaked (LSPs) blazars with 
the Vp eak < 10 14 ' 5 Hz [3j. Anton and Brown [4] also considered very low energy peaked BL Lacs (VLBL), 
with the synchrotron component peaking in IR/mm band Vp eak < 10 13 Hz. 

According to the works on synchrotron peak frequency estimation, only about 10% of blazars have i/* eafe < 
10 13 Hz, especially they are rare among BL Lacs. For example, in [5] only 3% of their sample could be 
classified as VLBL. According to the BLcat [6] (includes more than 300 BL Lacertae type objects, observed 
with the RATAN-600 radio telescope), about 9% are VLBL among BL Lacs. The typical value of the zV eafc 
for FSRQs (flat-spectrum radio quasars) is 10 131 ± 01 Hz, for BL Lacs i+ ea(fc value is one order of magnitude 
greater [7]. 

Nowadays many ground-based and space telescopes observing AGNs during different monitoring pro¬ 
grammes, also in order to build broadband SED. It is crucial to analyse quasi-simultaneous SED because of 
variable nature of blazars. The shape and positions of the components of SED could change during active 
state of the blazars, which will lead to misclassification. For example, synchrotron peak frequency shifted 
from 10 12 to 10 14 Hz during the flare in a gamma-ray band for object 4C+49.22 [8]. This object could be 
considered as a transition type, evolving from FSRQ to BL Lac [9, 10]. A less dramatic case of v* eak value 
change was detected in PKS 1510-089, that demonstrated shift from 1.5 x 10 13 to 6.5 x 10 13 Hz [11]. 

In this paper we will use quasi-simultaneous observations in optical (B, V, R bands) and radio (4.8-21.7 
GHz) domain to estimate z+ eafc values for six blazars, preliminarily defined by us as a candidates to the very 
low synchrotron peaked sources. Throughout this paper we will refer to our objects as very low synchrotron 
peaked blazars — VLSP. 
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2. THE SAMPLE OF VLSP BLAZARS 

The blazars monitoring list at the RATAN-600 radio telescope [12-14] includes different types of objects 
(FSRQs and BL Lacs), the BL Lacertae type objects list could be find in the on-line catalogue 1 [6]. We 
selected six blazars — VLSP candidates from our list to observe them with the Zeiss-1000 telescope at the 
same time with the RATAN-600 observations during year 2014. We made a preliminarily estimation of the 
Vpeak using tli e ASDC SED Builder Tool 2 . It is useful tool to create multi-waveband SED and made an 
polynomial approximation of the different regions of the SED to calculate the Vp eak value. In Table 1 we 
present: source name, right ascension (RA) and declination (DEC) (for the J2000.0 epoch), redshift value 3 , 
and preliminarily estimated i/® eofe value. 


3. OBSERVATIONS AND DATA REDUCTION 
3.1. Radio data 

We used the flux densities measured at five frequencies (2.3, 4.8, 7.7, 11.2 and 21.7 GHz) in 2014 February, 
March, July and December at the RATAN-600 radio telescope. Each source was observed from 5 to 12 times 
during this period 4 . The experimental data were processed with the modules of the FADPS (Flexible 
Astronomical Data Processing System) standard reduction package by [15]. The processing methods are 
described in the paper by [16]. The following 12 flux density calibrators (standard and RATAN’s traditional 
ones) were used to get the coefficients of antenna elevation: 3C48, 3C138, 3C147, 3C161, 3C286, 3C295, 
3C309.1, NGC7027, J0237—23, J1154—35, J1347+12 and J0410+76. Measurements of some calibrators were 
corrected for angular size and linear polarization, following the data summarized in [17] and [18], respectively. 
The detection limit for the RATAN-600 single sector is approximately 8 mJy (integration time is about 3 
s) under good weather conditions at the frequency of 4.8 GHz and at an average antenna elevation (S ~ 
42°). Standard errors in determining the flux density for these data are 3-5% at 4.8 u 7.7 GHz, 4-10% at 
11.2 GHz, and 10-14% at 21.7 GHz. All sources had relatively strong flux levels at radio frequencies with 
a signal-to-noise ratio S/N > 6. In Table 2, we list the radio flux densities at different frequencies used in 
this work. 


3.2. Optical data 

Photometric optical observations were performed in 2014 February, May, July and November at 1-meter 
Zeiss-1000 telescope, equipped with liquid nitrogen-cooled CCD camera EEV 42-40 (2048 x 2048 pixels). A 
standard B , V (Johnson) and R (Cousins) filters were used. Exposure times, depending upon the brightness 
of the source and weather conditions, were 30-300 s long 5 . The average size of the objects was 2-3"under 
good astroclimatic conditions. All the data reduction procedures were conducted using standard methods in 
the MaxlmDL. Reference stars were chosen with a brightness comparable to the target star, located in the 
same field of the CCD matrix. Thus, atmospheric extinction was the same for the reference and target star 
and did not require additional accounting. Correction for interstellar extinction was made using values from 
NED. The standard errors in determining the magnitude are 0.8% in B , 0.7% in V, and 0.4% in R filter. We 
used standard formula to convert magnitudes to flux densities: 


mi — m 2 = —2.5 log 


fi 

f2 


there mi, m 2 — are magnitudes of the target and reference stars, /i ,/2 — their flux densities, respectively. 
If we take the calibration flux density value (/o), at which magnitude is equal to zero (m = 0), i.e. write 


1 http://www.sao.ru/blcat/ 

2 http://tools.asdc.asi.it 

3 data from the NASA/IPAC (NED) http://ned.ipac.caltech.edu/ 

4 but for the [HB89] 1308+326 measurements available only at 7.7 GHz 

5 unfortunately the PKS 0528+134 was too faint in B and V filters to estimate its magnitude 
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down 


m 2 = 0, f 2 = jo, 


then 


and flux density will be 


mi = —2.5 log 


h 

/o’ 


— m -I 

fi = 10 2 5 jo- 

The zeropoint flux density values /o are taken from [19] and presented in Table 3 for each filter. The results 
of the photometry are presented in Table 4. 


4. RESULTS 

We plotted spectral energy distributions of our VLSP candidates using our measurements in radio and 
optical domain (flux densities presented in Table 2 and Table 4) in Fig 1. Most of the observed radio to 
optical (and in some cases X-ray) emission from blazars is due to synchrotron radiation of the relativistic 
electrons moving in a magnetic field [20-23]. Its spectrum described by a power low [24] and second or 
third degree polynomial is usually used to describe form of the synchrotron component of the spectrum (for 
example to find its maximum) [5, 25-28]. 

The synchrotron component of the SED (10 8 — 10 15 Hz) in our case was fitted with a parabolic function: 

log ( vF) = A + B{log v) + C(log u) 2 

We used the OriginLab (software for data analysis and graphing) to estimate synchrotron peak frequency. 
The results of z/® eofc estimation based on our almost concurrent and archival non-concurrent data for our 
sample are presented in Table 5. Comments on individual sources provided below. 

PKS 0336-01. This source has been classified as a flat-spectrum radio quasar (FSRQ) [29], with the 
dominating radio core. The optical light curve presented at the St. Petersburg State University (SPbSU) 
virtual observatory 6 shows decreasing of the R band magnitude at the end of the 2014: from ~15.7 to ~17.2. 
In October 2014 a GeV gamma-ray flare of this object was detected [30], but at the same time there was no 
activity in optical band [31]. Interestingly, after couple of weeks NIR flare of this blazar was detected [32]. 
According to the Owens Valley Radio Observatory (OVRO) monitoring 7 , source was stable at the radio band 
at the end of 2014 (flux density at 15 GHz ~2.4 Jy). 

Our observations were performed in November 19 and in December 7 12 in optical and radio band, 
respectively. We determined log y peak = 13.67 (four points, the reduced chi-square y 2 = 0.0092), which is 
somewhat higher than archival value log z/* eQfc = 12.68. 

PKS 0446+11. This is BL Lacertae type object [33], with almost featureless optical spectra [34], have a 
flat radio spectrum [29]. The source was in a bright state at radio band in early 2014, according to OVRO 
monitoring programme. Note that the radio light curve of this source characterized with a slow rise (and 
decline) during activity, maximum of a brightness is reached usually in a scale of about one year. The flux 
density variation is less than 10% in a scale of 1-2 weeks. 

Our observations were performed in February 2 in optical band and during February 12 March 8 in 
radio band. From the SED of this object (see Fig 1) it could be seen that the positions of the points in 
optical band influenced the parabolic fit, so that maximum became higher than it could be expected. We 
determined log v peak = 15.30 (six points, y 2 = 0.108), while estimated value from archival data is log v peak 
= 12.80. 

PKS 0528+134. This blazar has a radio jet [35] and classified as a FSRQ [29]. 

We observed this source in February 2 in optical band and during February 12 - March 8 in radio band. 
This source was generally very faint in optical band, and we obtained only three images in R filter (3 


6 http://lacerta.astro.spbu.ru/ 

' http://www.astro.caltech.edu/ovroblazars/data/data.php 
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expositions, 300 s each), we failed to accumulate sufficient amount of photons in B and V filters to calculate 
its magnitude. Later in that year, in 2014 November, we tried again but it was still in a faint state (less than 
19 mag in R filter). We estimated log i/® eafc = 13.77 (four points, y 2 = 0.0017); according to the archival 
data log v s peak = 11.96. 

[HB89] 1308+326. It is BL Lacertae type object, with weak emission lines, host galaxy is not de¬ 
tectable [36, 37]. A strong flare was detected in gamma-ray band in the 2014 April [38], but following optical 
observations did not reveal activity [39]. At the first half of 2014 this blazar was in active state in radio 
band (flux density at 15 GHz ~2 Jy, according to the OVRO monitoring programme). 

We observed this source in May 28 and five times in May 27-31 in optical and radio band, respectively. 
We estimated log i/® eofc = 12.79 (seven points, % 2 = 0.0025), which is in a good agreement with the value 
obtained using archival data — log Vp eak = 12.68. 

3C 345. This source has been classified as a FSRQ [29], with the dominating radio core [40]. It was in a 
quiescent state in optical band in 2014, slightly brighter than 18 mag in R filter (according to the SPbSU 
virtual observatory); slight fall in brightness in radio domain, which began in 2013, has continued — from 
~6 to ~5 Jy (OVRO data). 

Our observations were performed in July 24 in optical band and two weeks earlier, in July 7 12 in radio 
band. We estimated log ;/® eafc = 12.56 (four points, y 2 = 0.0079). The log i/* eofe = 12.78, according to the 
archival data. It is the second source in our sample for which our estimated value of the synchrotron peak 
is very close to that determined from archival non-concurrent data. 

PKS 2230+11. This source is known as a typical blazar [41] and as a quasar with the high optical 
polarization (>3%) [33]; a radio jet is also observed in the object [35]. The significant SED variation was 
detected during multi-waveband monitoring of its flare in 2005 [41]. One can approximately estimate the z/* eafc 
value decrease of about one order, with the spread of activity to the low-frequency region of the spectrum, 
from the SED presented in that work [41]. There was a brightening of about one order of magnitude in R 
band at the end of the 2014, according to the light curve presented in the SPbSU virtual observatory web 
page (in comparison with its stable state in the middle of 2014). It is possible that this activity is linked 
with the flare in a gamma-ray band, registered in late October of the same year [42]. 

We observed this source in November 19 in optical band and three weeks later, during December 7-12, in 
radio band. We determined log Vp eak = 13.64 (four points, y 2 = 0.0005); log i/| eofc = 12.86, archival data. 


5. CONCLUSIONS 

We have presented new observational results at radio and optical bands for six blazars, which were selected 
to be VLSP candidates (z 'p eak < 10 13 Hz) according to the non-concurrent archival data. We estimated 
synchrotron peak frequency Vp eak using quasi-simultaneous observations at the Zeiss-1000 and the RATAN- 
600 telescopes, the difference in observation epochs was not more than 2-3 weeks. 

We confirmed classification as VLSP for two sources: [HB89] 1308+326 h 3C345. Other four objects 
considered as candidates to VLSP (PKS 0336-01, PKS 0446+11, PKS 0528+134 and PKS 2230+11) have 
shown Vp eak > 10 13 Hz. We found synchrotron peak frequency values very close to the estimated from 
archival data for [HB89] 1308+326 and 3C345. 
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Figure 1 . Spectral energy distributions build using our RATAN-600 and Zeiss-1000 measurements for objects from our 

sample. 
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Table 1 . Sample of the VLSP candidates 


Name 

| RA (2000.0) 

DEC (2000.0) 

z 

log Vp ea k 

PKS 0336-01 

03 h 39 m 

30.9 s 

-01°46 36” 

0.852 

12.68 

PICS 0446+11 

04'M9 m 

07.7 s 

+11°21 , 29” 

1.207 

12.80 

PKS 0528+134 

05 h 30 m 

56.4 s 

+13°31 , 55” 

2.06 

11.96 

[HB89] 1308+326 

13 h 10 m 

28.6 s 

+32°20 , 44” 

0.998 

12.68 

3C 345 

16 ,l 42 m 

58.8 s 

+39°48 , 37” 

0.592 

12.78 

PKS 2230+11 

22 h 32 m 

36.4 s 

+ll°43 , 5l” 

1.037 

12.86 


Table 2. The flux density values at four frequencies measured with the RATAN-600 


Name 

Date 

F21.7GHZ, Jy 

Fll.2GHz, Jy 

Ft. 7 GHz, Jy 

Fa.sghz, Jy 

PKS 0336-01 
PKS 0446+11 
PKS 0528+134 
[HB89] 1308+326 
3C 345 

PKS 2230+11 

2014 Dec 7-12 
2014 Feb 12 - Mar 8 
2014 Feb 12 - Mar 8 
2014 May 27-31 
2014 July 7-12 
2014 Dec 7-12 

1.570 ± 0.194 
0.634 ± 0.078 
0.811 ± 0.116 

1.74 ± 0.07 
5.669 ± 0.538 

2.571 ± 0.315 

1.860 ± 0.180 
0.782 ± 0.031 
1.193 ± 0.048 
2.14 ± 0.11 
5.704 ± 0.285 
3.260 ± 0.241 

2.14 ± 0.11 

2.23 ± 0.096 
0.694 ± 0.021 
1.424 ± 0.043 
1.92 ± 0.23 
5.846 ± 0.234 
4.207 ± 0.181 


Table 3. Parameters of the optical filters and their calibration flux density values (/o for mag 


0) from [19] 


Band 

A, nm 

v, Hz 

fo, Jy 

B 

440 

6.81 x 10 14 

4260 

V 

550 

5.45 x 10 14 

3640 

R 

640 

4.68 x 10 14 

3080 


Table 4. Magnitudes at three bands, measured with the Zeiss-1000, and corresponding flux densities of the objects from our 

sample (this values before correction for interstellar extinction) 


Name 

Date 

B, mag 

V, mag 

R, mag 

B, mjy 

V, mjy 

R, mjy 

PKS 0336-01 
PKS 0446+11 
PKS 0528+134 
[HB89] 1308+326 
3C 345 

PKS 2230+11 

2014 Nov 19 
2014 Feb 2 
2014 Feb 2 
2014 May 28 
2014 July 24 
2014 Nov 19 

17.65 ± 0.12 
18.14 ± 0.15 

18.64 ± 0.09 
18.37 ± 0.26 
17.39 ± 0.10 

17.54 ± 0.09 
18.16 ± 0.13 

18.33 ± 0.12 
18.18 ± 0.20 
16.78 ± 0.07 

17.21 ± 0.05 
18.61 ± 0.10 
19.27 ± 0.08 
17.70 ± 0.05 
17.76 ± 0.07 
16.50 ± 0.06 

0.36 ± 0.0025 
0.229 ± 0.0019 

0.145 ± 0.0007 
0.185 ± 0.0026 
0.457 ± 0.0025 

0.364 ± 0.0018 
0.206 ± 0.0015 

0.176 ± 0.001 
0.202 ± 0.0023 
0.734 ± 0.0031 

0.384 ± 0.0012 
0.106 ± 0.0006 
0.0576 ± 0.0002 
0.245 ± 0.0007 
0.231 ± 0.0009 
0.739 ± 0.0026 


Table 5. Measurements of the z/* eafc parameter for the investigated objects from our sample (at column two, designated 
SAO), estimated values of the vp eak using SED Builder Tool also provided (at column three, designated archive). Objects, 

confirmed as VLSP, marked with + (at column four) 


Name 

log J'peafc, SAO 

log J'peafc, archive 

VLSP 

PKS 0336-01 

13.67 

12.68 

- 

PKS 0446+11 

15.30 

12.80 

- 

PKS 0528+134 

13.77 

11.96 

- 

[HB89] 1308+326 

12.79 

12.68 

+ 

3C 345 

12.55 

12.56 

+ 

PKS 2230+11 

13.64 

12.86 

- 



































